An estimate for the location of QCD critical end point 
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It is proposed that a study of the ratio of shear viscosity to entropy density ^ as a function of 
the baryon chemical potential /is, and temperature T, provides a dynamic probe for the critical 
end point (CEP) in hot and dense QCD matter. An initial estimate from an elliptic flow excitation 
function gives /i^^^ - 150 - 180 MeV and Tcep - 165 - 170 MeV for the location of the the CEP. 
These values place the CEP in the range for "immediate" validation at RHIC. 

PACS numbers: PACS 25.75.Ld 



The phase boundaries and the critical end point (CEP) 
are fundamental characteristics of hot and dense nuclear 
matter [1]. The study of heavy ion collisions has been 
proposed [ ] as an avenue to search for these essential 
characteristics of the Quantum Chromodynamics (QCD) 
phase diagram i.e. the plane of temperature vs baryon 
chemical potential (T, /i^). 

A recent resurgence of experimental interest in the 
CEP has been aided by strong experimental and theoret- 
ical evidence for a crossover transition to the quark gluon 
plasma (QGP) in heavy ion collisions at the Relativist ic 
Heavy Ion collider (RHIC) [3, 4, 5, 6, 7, 8, 9, 10, 11]. 
Such a crossover, constitutes a necessary requirement, 
albeit insufficient, for locating the CEP. 

Several attempts have been made to provide theoreti- 
cal guidance on where to localize a search for the CEP in 
the QCD phase diagram [7, 12, 13, 14, 15]. The resulting 
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predictions for the critical values of temperature Tcep and 
baryon chemical potential /i^^, which locates the CEP, 
have not converged and now span a broad range. 

Therefore, recent plans for the experimental verifica- 
tion of the CEP have centered on energy scans with an 
eye toward accessing the broadest possible range oi jiB 
and T values [16, 17, 18, 19]. Fig. 1 reenforces the value 
of such energy scans; it shows the chemical freezeout val- 
ues oi jiB (top panel) and T, as a function of beam colli- 
sion energy -^^snn^ extracted via chemical fits to particle 
ratios [20] obtained at several accelerator facilities. This 
unprecedented reach in /i^ and T values, clearly indi- 
cate that the combined results from energy scans at the 
Facility for Anti-proton and Ion Research (FAIR), the 
Super Proton Synchrotron (SPS) and RHIC, will allow 
access to the full range oi jiB and T values necessary for 
a comprehensive CEP search. 

At the CEP (or close to it) anomalies can occur in a 
wide variety of dynamic and static properties. Anoma- 
lies in dynamic properties reflect a change in quantities 
such as the transport coefficients and relaxation rates, 
multi-time correlation functions and the linear response 
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FIG. 1: (Color online) y/sNN dependence of /is (top panel) 
and T obtained from chemical fits [20] to particle yield ratios 
obtained at different accelerator facilities as indicated. The 
solid lines are fits to these data. 



to time-dependent perturbations. All of these depend 
on the equations of motion, and are not simply deter- 
mined by the equilibrium distribution of the particles 
at a given instant of time. By contrast, static proper- 
ties are solely determined by the single-time equilibrium 
distribution which include thermodynamic coefficients, 
single-time correlation functions, and the linear response 
to time-independent perturbations. 

Critical fluctuations are thought to be one of the more 
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important static signals for locating the CEP ["^T]. Con- 
sequently, extensive studies of particle fluctuations have 
been made over a broad range of beam collision energies. 
To date, no definitive observation indicative of the CEP, 
has been reported. A comprehensive search for the CEP 
via dynamic variables is still lacking. 

In this letter, we argue that it is possible to locate the 
CEP via study of the jub and T dependence of the ratio 
of viscosity to entropy density (^) [22, 23, 24] and give an 
estimate for /i^^ and Tcep (i.e the location of the CEP) 
using existing data. 

The rationale for using ^ as a probe for the CEP is 
two fold. First, we observe that this QCD critical end- 
point belongs to the universality class of the 3d Ising 
model, ie. the same universality class for a liquid-gas 
system; here, it is important to recall that all mem- 
bers of a given universality class have "identical" critical 
properties. Second, we observe that ^, for atomic and 
molecular substances, exhibits a minimum of compara- 
ble depth for different isobars passing in the vicinity of 
the liquid-gas critical end point [22, 23, 24, 25]. Fig. 2 
illustrates this for H2O. When an isobar passes through 
the critical end point, the minimum forms a cusp at the 
reduced temperature "^""^"""p = 0; when it passes above 
the critical end point (i.e a pressure P above the criti- 
cal pressure Peep), a less pronounced minimum is found 
at a value slightly above ^ ''''^ = 0. For an isobar 
passing just below the critical pressure, the minimum 
is found at ^^Jj""^ < (liquid side) but is accompanied 
by a discontinuous change across the phase transition. 
Thus, for a range of reduced temperatures, the average 
value (47r(^)) can be seen to grow rapidly for isobars 
passing through the critical end point and just below it. 
This is illustrated in the inset of Fig. 2 where (47r(^)) 
is plotted vs pressure for the reduced temperature range 
^^^^ = - 0.3. Fig. 2 shows that the CEP is signaled 

^ cep 

by a minimum at ~ ''''^ ~ 0, in the dependence of 47r(-) 

-L cep ^ S ' 

on the reduced temperature, as well as a sharp increase 
in (47r(f )) vs P for ^^^^ > 0. 

In analogy to the observations for atomic and molec- 
ular substances, one expects a range of trajectories, in 
the (T, /i^) plane for decaying nuclear systems, to show 
^ minima with a possible cusp at the critical end point 
(Tcep, Mb^)- That is, for /i^ = /i^^ the - minimum is ex- 
pected at the reduced temperature ^ ''''^ =0; for other 
values of fiB with associated critical temperature Tc not 
too far from Tcep, the dependence of 47r(^) on ^^Jj""^ is 
also expected to grow stronger as /i^ is increased from 
an initially small value up to /i^ ^ Mb^- Iiideed, re- 
cent calculations for different types of phase transitions 
(first-order, second-order and a crossover) suggest a rapid 
change in the value of 7^/5 in the vicinity of the CEP [26]. 

For a given value of /i^ , a hot nuclear system for which 
^"^''''^ > will sample the full range of ^ values to 

-L cep S 



give an average, as it evolves toward the | minimum. 
Consequently, one expects an increase of (47r(|)) with 
increasing /i^, punctuated by a relatively rapid increase 
for jiiB values slightly above /i^^. The latter would be 
comparable to the rapid increase in (47r(^)) observed for 
II2O (inset in Fig. 2) when P is lowered a little below the 
critical pressure. 

Therefore, the extraction of (47r(^)) as a function of T 
and /i B from experimental data, can serve as a constraint 
for the location of the CEP. Such extractions are possible 
from an elliptic fiow excitation function measurement be- 
cause a sizable change in (47r(^)) is expected to lead to a 
measurable suppression of the magnitude of elliptic fiow. 
It could even serve to invalidate the currently observed 
universal scaling patterns [27, 28, 29]. 

In recent work [23], we have used elliptic fiow mea- 
surements to obtain the estimates (47r(^)) ~ 1.3 [30] and 
(T) ~ 165 MeV for hot and dense matter [31] produced in 
Au+Au collisions {^/s^ = 200 GeV or /i^ ^ 24 MeV) 
at RHIC. A comparison of this ^ value to those calcu- 
lated for a meson-gas for T < Tf [24], and the QGP for 
T > Tf (Tf - 170 MeV [ ]), gave a good indication for 
the expected minimum (for T close to Tcep) in the plot 

of 47r(-) vs ^"q"^ . We therefore use this observation 
as a basis for the estimate Tcep ~ 165 — 170 MeV. This 
estimate is similar to the chemical freeze-out tempera- 
ture for a broad range of collision energies (see botom 
panel of Fig. 1). This constancy of the freeze-out tem- 
perature (T ~ 165 MeV) may be a further indication 
that chemical freeze-out occurs at, or close to Tcep for 
V^i^ - 17 - 200 GeV. 

The value (47r(^)) ~ 1.3, achieved in Au+Au collisions 
at s/snn = 200 GeV, is rather close to the conjectured 
lower bound of 47r(^) = 1.0. Consequently, one can con- 
clude that, ioi jiB ^ 24 MeV, the hot expanding system 
spends a considerable portion of its dynamics in the re- 
gion of low ^; an optimal situation being, the system 
stays at low | and then very quickly freezes out at or 
close to Tcep. Such a trajectory would be tantamount to 
the isobaric trajectory above the critical pressure, shown 
for H2O in Fig. 2 (triangles). For other trajectories with 
11 B close to, or slightly above /i^^, significant collective 
motion is expected to develop as the system evolves to- 
ward freeze-out with higher ^ values. Consequently, the 
dependence of elliptic fiow on jiB is of interest. 

Figure 3 shows a differential elliptic fiow {V2) excita- 
tion function for charged hadrons, {pr) = 0.65 GeV/c, 
measured in 13 — 26% central Au-hAu and Pb-hPb col- 
lisions [33]. For the range of collision energies -^Z ^ 
17 — 200 GeV, Fig. 1 indicates an essentially constant 
freeze-out temperature T ~ 165 MeV for the values 
fiB ~ 25 — 250 MeV. Therefore, these V2 measurements 
(for -^snn ~ 17 — 200 GeV) result from excited systems 
which all evolve toward our assumed value of Tcep, albeit 
with different fiB values. 
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2: (Color online) 47r(^) vs the reduced temperature 
££p for H2O. Results are shown for an isobar at the critical 

cep 

pressure (Peep = 22.6 MPa) and one above (below) it, as indi- 
cated. The inset shows (47r(^)) for the reduced temperature 
range 0.0 — 0.3. The data are taken from Ref. [22]. 



Figure 3 shows that V2 is essentially constant for 
y/sNN ^ 62 — 200 GeV. The energy density is estimated 
to decrease by ~ 30% as the beam collison energy is 
reduced from ^snn ^ 200 GeV to -^snn ^ 62 GeV. 
Therefore we interpret this constancy of V2 as an indi- 
cation that (i) the equation of state associated with the 
crossover transition to the QGP is soft, and (ii) that 
is relatively small for the /i b values corresponding to this 
collision energy range. That is, these jub values are sig- 
nificantly smaller than /i^^. 

For -^snn ~ 18 GeV Fig. 3 shows that V2 decreases 
by almost 50%, compared to the value for ^/sNN ^ 
62 — 200 GeV. Here, it is important to point out that 
the mean transverse energy per particle is essentially the 
same for collision energy range ^snn ~ 17 — 200 GeV 
and the estimated Bjorken energy density is ~ 5.4 and 
3.2 GeV/fm^ for = 200 GeV [ ] and ^^s^ = 

17 GeV [ ] respectively, i.e the energy density change 
from -^snn ^ 62 GeV to -^snn ~ 17 GeV is not very 
large. Thus, the initial temperature of the high energy 
density matter created in collisions at ^snn ^ 62 GeV 
and ^snn ~ 17 GeV are not drastically different, and 
the fraction of the elliptic flow generated during the (dis- 
sipative) hadronic phase [ou] is expected to be qualita- 
tively similar. 

A reduction in collision energy from -^snn 62 GeV 
to -s/snn ~ 17 GeV leads to a significant increase (more 
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FIG. 3: (Color online) Flow excitation function. The data 
is obtained from Ref. [33]. The inset shows the schematic 
variation of (47r(^)) vs /is. 



than a factor of two) in the value of /i^. Recent cal- 
culations .] also indicate that, in the hadronic phase, 
^ decreases with increasing /i^. Therefore, a significant 
part of the reduction in V2 observed as the collision energy 
is reduced from ^/snn ^ 62 GeV to ^snn ~ 18 GeV, 
could be a manifestation of the expected increase in ( ^ ) 
for values oi jis^ Mb^- 

To estimate /i^^, we assume a smooth transition in the 
magnitude of V2 over the range 62 ^ ^Jsnn ^18 GeV, 
which is not yet measured. The inset in Fig. 3 gives a 
schematic illustration of the expected change of (47r(^)) 
with (iB over this collision energy range. We use the 
"knee" in the extrapolated values for V2 over the range 
62 ^ y/SNN ^ 18 GeV to obtain the estimate /i^^ ~ 
150 — 180 MeV. A similar estimate was obtained by evalu- 
ating the /i B dependence of ( |) following the procedures 
outlined in Refs. [22, 23, 38, 39], followed by interpola- 
tion to unmeasured iib values. 

Given these values of Tcep and /i^^, we expect the 
value ^47r(^)), extracted from flow measurements per- 
formed at -^snn ^ 40 and 30 GeV, to be significantly 
larger than those obtained from similar measurements 
performed over the range ^snn ~ 62 — 200 GeV. This 
change should also be reflected in the onset of a de- 
crease of ^2, a possible violation of the universal scal- 
ing patterns [27, 28, 29] observed for measurements at 
■^snn ~ 62 — 200 GeV and a measurable increase in V2 
fluctuations. 

In summary we have argued that experimental as- 
sessment of (1) as a function of fiB and T provides a 
good dynamic observable for constraining the critical end 
point of hot QCD matter. A first estimate for the CEP 
from flow data indicate the values Tcep ~ 165 — 170 and 
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150 — 180 MeV. Interestingly, our estimate is 
in good agreement with the prediction of Gavai et al. 
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[14], obtained from lattice QCD simulations with real- 
istic (about 1.7 times) pion masses and large volumes. 
This estimate also places the CEP in the range for direct 
validation at RHIC via an energy scan. An initial mea- 
surement at ^snn ^ 40 and 30 GeV would give sufficient 
information on where to focus more detailed attention. 
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